(* Content-type: application/vnd.wolfram.mathematica *)

(*** Wolfram Notebook File ***)
(* http://www.wolfram.com/nb *)

(* CreatedBy='Mathematica 8.0' *)

(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[       157,          7]
NotebookDataLength[     26229,        820]
NotebookOptionsPosition[     23773,        737]
NotebookOutlinePosition[     24110,        752]
CellTagsIndexPosition[     24067,        749]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{
Cell[BoxData[
 RowBox[{
  RowBox[{"M", "[", "x_", "]"}], " ", ":=", " ", 
  RowBox[{
   RowBox[{
    RowBox[{"a", " ", "/", "3"}], "*", 
    RowBox[{"(", 
     RowBox[{
      RowBox[{
       RowBox[{"x", "^", "4"}], "/", "4"}], " ", "-", " ", 
      RowBox[{
       RowBox[{"L", "^", "3"}], "*", "x"}], " ", "+", " ", 
      RowBox[{"3", "*", 
       RowBox[{
        RowBox[{"L", "^", "4"}], "/", "4"}]}]}], ")"}]}], " ", "+", " ", 
   RowBox[{
    RowBox[{"b", "/", "2"}], "*", 
    RowBox[{"(", 
     RowBox[{
      RowBox[{
       RowBox[{"x", "^", "3"}], "/", "3"}], " ", "-", " ", 
      RowBox[{
       RowBox[{"L", "^", "2"}], "*", "x"}], "+", 
      RowBox[{"2", "*", 
       RowBox[{
        RowBox[{"L", "^", "3"}], "/", "3"}]}]}], ")"}]}], " ", "+", " ", 
   RowBox[{"c", "*", 
    RowBox[{"(", 
     RowBox[{
      RowBox[{
       RowBox[{"x", "^", "2"}], "/", "2"}], " ", "-", " ", 
      RowBox[{"L", "*", "x"}], " ", "+", " ", 
      RowBox[{
       RowBox[{"L", "^", "2"}], "/", "2"}]}], ")"}]}], " ", "+", " ", 
   RowBox[{"F", 
    RowBox[{"(", 
     RowBox[{"L", "-", "x"}], ")"}], " "}]}]}]], "Input",
 CellChangeTimes->{{3.5363317274022408`*^9, 3.5363318843186903`*^9}, {
  3.536332021384294*^9, 3.5363320503540983`*^9}, {3.5365820062182493`*^9, 
  3.5365820092667294`*^9}, {3.5371131297091513`*^9, 3.537113131984202*^9}, {
  3.537113211665009*^9, 3.537113225893674*^9}, {3.537113308308111*^9, 
  3.537113345968528*^9}, {3.5371139148274097`*^9, 3.537113916615883*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"M", "[", "x", "]"}]], "Input",
 CellChangeTimes->{{3.536331887424693*^9, 3.536331895546067*^9}, {
  3.53711265107169*^9, 3.5371126520345697`*^9}, {3.537113120661242*^9, 
  3.537113122156863*^9}, {3.5371133115501738`*^9, 3.537113312799591*^9}}],

Cell[BoxData[
 RowBox[{"F", " ", 
  RowBox[{"(", 
   RowBox[{"L", "-", "x"}], ")"}]}]], "Output",
 CellChangeTimes->{
  3.5363318966256647`*^9, 3.5363320671111603`*^9, 3.536335457846232*^9, 
   3.537106187307394*^9, 3.537111302926835*^9, 3.537112586914768*^9, 
   3.537112652502781*^9, {3.5371131225579967`*^9, 3.537113136261718*^9}, 
   3.537113283284693*^9, 3.537113313558086*^9, 3.537113348505249*^9, 
   3.537113919173917*^9, 3.5371142447642727`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Integrate", "[", 
  RowBox[{
   RowBox[{
    RowBox[{"M", "[", "x", "]"}], "*", 
    RowBox[{"M", "[", "x", "]"}]}], ",", 
   RowBox[{"{", 
    RowBox[{"x", ",", "xi", ",", "xj"}], "}"}]}], "]"}]], "Input",
 CellChangeTimes->{{3.536331974519031*^9, 3.536331995320217*^9}, {
  3.5363320929018373`*^9, 3.5363320950633783`*^9}, {3.53633377210667*^9, 
  3.5363337758432426`*^9}}],

Cell["\<\
(-(1/3))*F^2*(xi - xj)*(3*L^2 + xi^2 + xi*xj + xj^2 - 3*L*(xi + xj))\
\>", "Output",
 GeneratedCell->False,
 CellAutoOverwrite->False,
 CellChangeTimes->{
  3.536331997395433*^9, {3.536332078849463*^9, 3.5363320985565357`*^9}, 
   3.536333780427396*^9, 3.536333908994501*^9, {3.536334635964037*^9, 
   3.536334652621677*^9}, 3.536335459214384*^9, 3.537106192695814*^9, 
   3.537111308259521*^9, 3.5371125895770397`*^9, {3.537113445005517*^9, 
   3.537113460369835*^9}, {3.5371139232438498`*^9, 3.53711393593711*^9}, 
   3.537114246457909*^9, 3.5371143106516247`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"M", "[", "x", "]"}]], "Input",
 CellChangeTimes->{{3.536333370975418*^9, 3.536333373873847*^9}}],

Cell[BoxData[
 RowBox[{"F", " ", 
  RowBox[{"(", 
   RowBox[{"L", "-", "x"}], ")"}]}]], "Output",
 CellChangeTimes->{3.5363333746719933`*^9, 3.5363354593285437`*^9, 
  3.537106192801059*^9, 3.537111308433769*^9, 3.537112589668961*^9, 
  3.537113578749679*^9, 3.537113976045561*^9, 3.5371142573682547`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"M", "'"}], "[", "x", "]"}]], "Input",
 CellChangeTimes->{{3.536333375952282*^9, 3.536333377753984*^9}}],

Cell[BoxData[
 RowBox[{"-", "F"}]], "Output",
 CellChangeTimes->{3.536333378326026*^9, 3.53633545943471*^9, 
  3.537106192874975*^9, 3.5371113084935503`*^9, 3.537112589753985*^9, 
  3.5371135801026077`*^9, 3.5371139767730227`*^9, 3.537114257850852*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"M", "''"}], "[", "x", "]"}]], "Input",
 CellChangeTimes->{{3.536332334703062*^9, 3.5363323609830933`*^9}}],

Cell[BoxData["0"], "Output",
 CellChangeTimes->{3.536332362096579*^9, 3.5363354595110607`*^9, 
  3.53710619295928*^9, 3.5371113085649643`*^9, 3.537112589839658*^9, 
  3.537113581191889*^9, 3.5371139776525583`*^9, 3.537114258576028*^9}]
}, Open  ]],

Cell[BoxData[
 RowBox[{"\[IndentingNewLine]", 
  RowBox[{
   RowBox[{
    RowBox[{"H1", "[", "x_", "]"}], " ", ":=", " ", 
    RowBox[{"1", " ", "-", " ", 
     RowBox[{"3", "*", 
      RowBox[{
       RowBox[{
        RowBox[{"(", 
         RowBox[{"x", "-", "xi"}], ")"}], "^", "2"}], "/", 
       RowBox[{
        RowBox[{"(", 
         RowBox[{"xj", "-", "xi"}], ")"}], "^", "2"}]}]}], " ", "+", " ", 
     RowBox[{"2", "*", 
      RowBox[{
       RowBox[{
        RowBox[{"(", 
         RowBox[{"x", "-", "xi"}], ")"}], "^", "3"}], "/", 
       RowBox[{
        RowBox[{"(", 
         RowBox[{"xj", "-", "xi"}], ")"}], "^", "3"}]}]}]}]}], " ", 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"H2", "[", "x_", "]"}], " ", ":=", " ", 
    RowBox[{
     RowBox[{"(", 
      RowBox[{"x", "-", "xi"}], ")"}], " ", "-", " ", 
     RowBox[{"2", "*", 
      RowBox[{
       RowBox[{
        RowBox[{"(", 
         RowBox[{"x", "-", "xi"}], ")"}], "^", "2"}], "/", 
       RowBox[{"(", 
        RowBox[{"xj", "-", "xi"}], ")"}]}]}], " ", "+", " ", 
     RowBox[{
      RowBox[{
       RowBox[{"(", 
        RowBox[{"x", "-", "xi"}], ")"}], "^", "3"}], "/", 
      RowBox[{
       RowBox[{"(", 
        RowBox[{"xj", "-", "xi"}], ")"}], "^", "2"}]}]}]}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"H3", "[", "x_", "]"}], " ", ":=", " ", 
    RowBox[{
     RowBox[{"3", "*", 
      RowBox[{
       RowBox[{
        RowBox[{"(", 
         RowBox[{"x", "-", "xi"}], ")"}], "^", "2"}], "/", 
       RowBox[{
        RowBox[{"(", 
         RowBox[{"xj", "-", "xi"}], ")"}], "^", "2"}]}]}], " ", "-", " ", 
     RowBox[{"2", "*", 
      RowBox[{
       RowBox[{
        RowBox[{"(", 
         RowBox[{"x", "-", "xi"}], ")"}], "^", "3"}], "/", 
       RowBox[{
        RowBox[{"(", 
         RowBox[{"xj", "-", "xi"}], ")"}], "^", "3"}]}]}]}]}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"H4", "[", "x_", "]"}], " ", ":=", " ", 
    RowBox[{
     RowBox[{
      RowBox[{"-", 
       RowBox[{
        RowBox[{"(", 
         RowBox[{"x", "-", "xi"}], ")"}], "^", "2"}]}], "/", 
      RowBox[{"(", 
       RowBox[{"xj", "-", "xi"}], ")"}]}], " ", "+", " ", 
     RowBox[{
      RowBox[{
       RowBox[{"(", 
        RowBox[{"x", "-", "xi"}], ")"}], "^", "3"}], "/", 
      RowBox[{
       RowBox[{"(", 
        RowBox[{"xj", "-", "xi"}], ")"}], "^", "2"}]}]}]}]}]}]], "Input",
 CellChangeTimes->{{3.536332543324222*^9, 3.536332646719583*^9}, {
  3.5363326947001257`*^9, 3.5363327624181757`*^9}, {3.536333129515232*^9, 
  3.5363331964174747`*^9}, {3.536337456507444*^9, 3.5363374587287292`*^9}, {
  3.536337562251939*^9, 3.5363375667725277`*^9}}],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{"H1", "[", "x", "]"}], "\[IndentingNewLine]", 
 RowBox[{"H2", "[", "x", "]"}], "\[IndentingNewLine]", 
 RowBox[{"H3", "[", "x", "]"}], "\[IndentingNewLine]", 
 RowBox[{"H4", "[", "x", "]"}]}], "Input",
 CellChangeTimes->{{3.536337461899653*^9, 3.536337483102229*^9}}],

Cell[BoxData[
 RowBox[{"1", "+", 
  FractionBox[
   RowBox[{"2", " ", 
    SuperscriptBox[
     RowBox[{"(", 
      RowBox[{"x", "-", "xi"}], ")"}], "3"]}], 
   SuperscriptBox[
    RowBox[{"(", 
     RowBox[{
      RowBox[{"-", "xi"}], "+", "xj"}], ")"}], "3"]], "-", 
  FractionBox[
   RowBox[{"3", " ", 
    SuperscriptBox[
     RowBox[{"(", 
      RowBox[{"x", "-", "xi"}], ")"}], "2"]}], 
   SuperscriptBox[
    RowBox[{"(", 
     RowBox[{
      RowBox[{"-", "xi"}], "+", "xj"}], ")"}], "2"]]}]], "Output",
 CellChangeTimes->{{3.536337479909708*^9, 3.536337483740776*^9}, 
   3.536337570295877*^9, 3.5371061930982857`*^9, 3.537111308742208*^9, 
   3.5371125899600563`*^9, 3.537113986618389*^9, 3.537114260473186*^9}],

Cell[BoxData[
 RowBox[{"x", "-", "xi", "+", 
  FractionBox[
   SuperscriptBox[
    RowBox[{"(", 
     RowBox[{"x", "-", "xi"}], ")"}], "3"], 
   SuperscriptBox[
    RowBox[{"(", 
     RowBox[{
      RowBox[{"-", "xi"}], "+", "xj"}], ")"}], "2"]], "-", 
  FractionBox[
   RowBox[{"2", " ", 
    SuperscriptBox[
     RowBox[{"(", 
      RowBox[{"x", "-", "xi"}], ")"}], "2"]}], 
   RowBox[{
    RowBox[{"-", "xi"}], "+", "xj"}]]}]], "Output",
 CellChangeTimes->{{3.536337479909708*^9, 3.536337483740776*^9}, 
   3.536337570295877*^9, 3.5371061930982857`*^9, 3.537111308742208*^9, 
   3.5371125899600563`*^9, 3.537113986618389*^9, 3.537114260475786*^9}],

Cell[BoxData[
 RowBox[{
  RowBox[{"-", 
   FractionBox[
    RowBox[{"2", " ", 
     SuperscriptBox[
      RowBox[{"(", 
       RowBox[{"x", "-", "xi"}], ")"}], "3"]}], 
    SuperscriptBox[
     RowBox[{"(", 
      RowBox[{
       RowBox[{"-", "xi"}], "+", "xj"}], ")"}], "3"]]}], "+", 
  FractionBox[
   RowBox[{"3", " ", 
    SuperscriptBox[
     RowBox[{"(", 
      RowBox[{"x", "-", "xi"}], ")"}], "2"]}], 
   SuperscriptBox[
    RowBox[{"(", 
     RowBox[{
      RowBox[{"-", "xi"}], "+", "xj"}], ")"}], "2"]]}]], "Output",
 CellChangeTimes->{{3.536337479909708*^9, 3.536337483740776*^9}, 
   3.536337570295877*^9, 3.5371061930982857`*^9, 3.537111308742208*^9, 
   3.5371125899600563`*^9, 3.537113986618389*^9, 3.537114260478257*^9}],

Cell[BoxData[
 RowBox[{
  FractionBox[
   SuperscriptBox[
    RowBox[{"(", 
     RowBox[{"x", "-", "xi"}], ")"}], "3"], 
   SuperscriptBox[
    RowBox[{"(", 
     RowBox[{
      RowBox[{"-", "xi"}], "+", "xj"}], ")"}], "2"]], "-", 
  FractionBox[
   SuperscriptBox[
    RowBox[{"(", 
     RowBox[{"x", "-", "xi"}], ")"}], "2"], 
   RowBox[{
    RowBox[{"-", "xi"}], "+", "xj"}]]}]], "Output",
 CellChangeTimes->{{3.536337479909708*^9, 3.536337483740776*^9}, 
   3.536337570295877*^9, 3.5371061930982857`*^9, 3.537111308742208*^9, 
   3.5371125899600563`*^9, 3.537113986618389*^9, 3.537114260480587*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{"H1", "''"}], "[", "x", "]"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"H2", "''"}], "[", "x", "]"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"H3", "''"}], "[", "x", "]"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"H4", "''"}], "[", "x", "]"}]}], "Input",
 CellChangeTimes->{{3.536332774578761*^9, 3.536332793960774*^9}, {
  3.5363348314134293`*^9, 3.53633483375905*^9}}],

Cell[BoxData[
 RowBox[{
  FractionBox[
   RowBox[{"12", " ", 
    RowBox[{"(", 
     RowBox[{"x", "-", "xi"}], ")"}]}], 
   SuperscriptBox[
    RowBox[{"(", 
     RowBox[{
      RowBox[{"-", "xi"}], "+", "xj"}], ")"}], "3"]], "-", 
  FractionBox["6", 
   SuperscriptBox[
    RowBox[{"(", 
     RowBox[{
      RowBox[{"-", "xi"}], "+", "xj"}], ")"}], "2"]]}]], "Output",
 CellChangeTimes->{3.5363331997739887`*^9, 3.536334836229802*^9, 
  3.536335460097742*^9, 3.536337573820105*^9, 3.537106193185505*^9, 
  3.537111308849001*^9, 3.53711259270304*^9, 3.5371135969117813`*^9, 
  3.537113988252898*^9, 3.5371142617877398`*^9}],

Cell[BoxData[
 RowBox[{
  FractionBox[
   RowBox[{"6", " ", 
    RowBox[{"(", 
     RowBox[{"x", "-", "xi"}], ")"}]}], 
   SuperscriptBox[
    RowBox[{"(", 
     RowBox[{
      RowBox[{"-", "xi"}], "+", "xj"}], ")"}], "2"]], "-", 
  FractionBox["4", 
   RowBox[{
    RowBox[{"-", "xi"}], "+", "xj"}]]}]], "Output",
 CellChangeTimes->{3.5363331997739887`*^9, 3.536334836229802*^9, 
  3.536335460097742*^9, 3.536337573820105*^9, 3.537106193185505*^9, 
  3.537111308849001*^9, 3.53711259270304*^9, 3.5371135969117813`*^9, 
  3.537113988252898*^9, 3.537114261790098*^9}],

Cell[BoxData[
 RowBox[{
  RowBox[{"-", 
   FractionBox[
    RowBox[{"12", " ", 
     RowBox[{"(", 
      RowBox[{"x", "-", "xi"}], ")"}]}], 
    SuperscriptBox[
     RowBox[{"(", 
      RowBox[{
       RowBox[{"-", "xi"}], "+", "xj"}], ")"}], "3"]]}], "+", 
  FractionBox["6", 
   SuperscriptBox[
    RowBox[{"(", 
     RowBox[{
      RowBox[{"-", "xi"}], "+", "xj"}], ")"}], "2"]]}]], "Output",
 CellChangeTimes->{3.5363331997739887`*^9, 3.536334836229802*^9, 
  3.536335460097742*^9, 3.536337573820105*^9, 3.537106193185505*^9, 
  3.537111308849001*^9, 3.53711259270304*^9, 3.5371135969117813`*^9, 
  3.537113988252898*^9, 3.537114261792254*^9}],

Cell[BoxData[
 RowBox[{
  FractionBox[
   RowBox[{"6", " ", 
    RowBox[{"(", 
     RowBox[{"x", "-", "xi"}], ")"}]}], 
   SuperscriptBox[
    RowBox[{"(", 
     RowBox[{
      RowBox[{"-", "xi"}], "+", "xj"}], ")"}], "2"]], "-", 
  FractionBox["2", 
   RowBox[{
    RowBox[{"-", "xi"}], "+", "xj"}]]}]], "Output",
 CellChangeTimes->{3.5363331997739887`*^9, 3.536334836229802*^9, 
  3.536335460097742*^9, 3.536337573820105*^9, 3.537106193185505*^9, 
  3.537111308849001*^9, 3.53711259270304*^9, 3.5371135969117813`*^9, 
  3.537113988252898*^9, 3.537114261794403*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"q", "[", "x", "]"}]], "Input",
 CellChangeTimes->{{3.536333061122662*^9, 3.536333062898327*^9}}],

Cell[BoxData[
 RowBox[{"q", "[", "x", "]"}]], "Output",
 CellChangeTimes->{3.5363330635656843`*^9, 3.536335460493999*^9, 
  3.536337574861721*^9, 3.5371061977673817`*^9, 3.53711130898847*^9, 
  3.537112594120775*^9, 3.537113600831686*^9, 3.53711399026128*^9, 
  3.5371142625992126`*^9}]
}, Open  ]],

Cell[BoxData[{
 RowBox[{
  RowBox[{"uM21", "[", "x_", "]"}], " ", ":=", " ", 
  RowBox[{
   RowBox[{"M", "[", "x", "]"}], "*", 
   RowBox[{
    RowBox[{"H1", "''"}], "[", "x", "]"}]}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"uM22", "[", "x_", "]"}], " ", ":=", " ", 
  RowBox[{
   RowBox[{"M", "[", "x", "]"}], "*", 
   RowBox[{
    RowBox[{"H2", "''"}], "[", "x", "]"}]}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"uM23", "[", "x_", "]"}], " ", ":=", " ", 
  RowBox[{
   RowBox[{"M", "[", "x", "]"}], "*", 
   RowBox[{
    RowBox[{"H3", "''"}], "[", "x", "]"}]}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"uM24", "[", "x_", "]"}], " ", ":=", " ", 
  RowBox[{
   RowBox[{"M", "[", "x", "]"}], "*", 
   RowBox[{
    RowBox[{"H4", "''"}], "[", "x", "]"}]}]}]}], "Input",
 CellChangeTimes->{{3.536332486562553*^9, 3.5363324950712357`*^9}, 
   3.53633253790725*^9, {3.536333220388173*^9, 3.536333249659178*^9}, {
   3.536334864560199*^9, 3.5363348667771187`*^9}, {3.536336497612015*^9, 
   3.536336502517456*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"\[IndentingNewLine]", 
  RowBox[{
   RowBox[{"V1", " ", "=", " ", 
    RowBox[{"Integrate", "[", 
     RowBox[{
      RowBox[{"uM21", "[", "x", "]"}], ",", 
      RowBox[{"{", 
       RowBox[{"x", ",", "xi", ",", "xj"}], "}"}]}], "]"}]}], 
   "\[IndentingNewLine]", 
   RowBox[{"V2", " ", "=", " ", 
    RowBox[{"Integrate", "[", 
     RowBox[{
      RowBox[{"uM22", "[", "x", "]"}], ",", 
      RowBox[{"{", 
       RowBox[{"x", ",", "xi", ",", "xj"}], "}"}]}], "]"}]}], 
   "\[IndentingNewLine]", 
   RowBox[{"V3", " ", "=", " ", 
    RowBox[{"Integrate", "[", 
     RowBox[{
      RowBox[{"uM23", "[", "x", "]"}], ",", 
      RowBox[{"{", 
       RowBox[{"x", ",", "xi", ",", "xj"}], "}"}]}], "]"}]}], 
   "\[IndentingNewLine]", 
   RowBox[{"V4", " ", "=", " ", 
    RowBox[{"Integrate", "[", 
     RowBox[{
      RowBox[{"uM24", "[", "x", "]"}], ",", 
      RowBox[{"{", 
       RowBox[{"x", ",", "xi", ",", "xj"}], "}"}]}], "]"}]}], 
   "\[IndentingNewLine]"}]}]], "Input",
 CellChangeTimes->{{3.536333269883795*^9, 3.536333332422559*^9}, 
   3.536335394909655*^9}],

Cell[BoxData[
 RowBox[{"-", "F"}]], "Output",
 CellChangeTimes->{{3.536335497499475*^9, 3.536335501058962*^9}, {
   3.5363365057059193`*^9, 3.536336530047544*^9}, 3.536337576650955*^9, 
   3.537106199584158*^9, 3.5371113105531273`*^9, 3.537112597144576*^9, {
   3.537113616142518*^9, 3.537113632867104*^9}, {3.537113993095098*^9, 
   3.5371140045146227`*^9}, 3.537114264011161*^9}],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"-", "F"}], " ", "L"}], "+", 
  RowBox[{"F", " ", "xi"}]}]], "Output",
 CellChangeTimes->{{3.536335497499475*^9, 3.536335501058962*^9}, {
   3.5363365057059193`*^9, 3.536336530047544*^9}, 3.536337576650955*^9, 
   3.537106199584158*^9, 3.5371113105531273`*^9, 3.537112597144576*^9, {
   3.537113616142518*^9, 3.537113632867104*^9}, {3.537113993095098*^9, 
   3.5371140045146227`*^9}, 3.537114264081296*^9}],

Cell[BoxData["F"], "Output",
 CellChangeTimes->{{3.536335497499475*^9, 3.536335501058962*^9}, {
   3.5363365057059193`*^9, 3.536336530047544*^9}, 3.536337576650955*^9, 
   3.537106199584158*^9, 3.5371113105531273`*^9, 3.537112597144576*^9, {
   3.537113616142518*^9, 3.537113632867104*^9}, {3.537113993095098*^9, 
   3.5371140045146227`*^9}, 3.537114264128085*^9}],

Cell[BoxData[
 RowBox[{
  RowBox[{"F", " ", "L"}], "-", 
  RowBox[{"F", " ", "xj"}]}]], "Output",
 CellChangeTimes->{{3.536335497499475*^9, 3.536335501058962*^9}, {
   3.5363365057059193`*^9, 3.536336530047544*^9}, 3.536337576650955*^9, 
   3.537106199584158*^9, 3.5371113105531273`*^9, 3.537112597144576*^9, {
   3.537113616142518*^9, 3.537113632867104*^9}, {3.537113993095098*^9, 
   3.5371140045146227`*^9}, 3.537114264188579*^9}]
}, Open  ]],

Cell[BoxData[{
 RowBox[{
  RowBox[{"phi1", "[", "x_", " ", "]"}], " ", ":=", " ", 
  RowBox[{"1", "-", 
   RowBox[{"x", "/", "l"}]}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"phi2", "[", "x_", "]"}], " ", ":=", " ", 
  RowBox[{"x", "/", "l"}]}]}], "Input",
 CellChangeTimes->{{3.5371062189459143`*^9, 3.537106248043092*^9}, {
  3.537111376787245*^9, 3.537111379192234*^9}}],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{"phi1", "[", "x", "]"}], "\[IndentingNewLine]", 
 RowBox[{"phi2", "[", "x", "]"}]}], "Input",
 CellChangeTimes->{{3.537106252234582*^9, 3.537106265032083*^9}}],

Cell[BoxData[
 RowBox[{"1", "-", 
  FractionBox["x", "l"]}]], "Output",
 CellChangeTimes->{3.5371062660218773`*^9, 3.5371113128675327`*^9, 
  3.5371113826123457`*^9, 3.537112600732087*^9}],

Cell[BoxData[
 FractionBox["x", "l"]], "Output",
 CellChangeTimes->{3.5371062660218773`*^9, 3.5371113128675327`*^9, 
  3.5371113826123457`*^9, 3.537112600734066*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{"xi", " ", "=", " ", "0"}], "\[IndentingNewLine]", 
 RowBox[{"xj", " ", "=", " ", "l", " "}]}], "Input",
 CellChangeTimes->{{3.53711140820868*^9, 3.537111414476397*^9}}],

Cell[BoxData["0"], "Output",
 CellChangeTimes->{3.5371114154271727`*^9}],

Cell[BoxData["l"], "Output",
 CellChangeTimes->{3.537111415428651*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{"C11", " ", "=", " ", 
  RowBox[{"Integrate", "[", 
   RowBox[{
    RowBox[{
     RowBox[{"phi1", "[", "x", "]"}], "*", 
     RowBox[{
      RowBox[{"H1", "''"}], "[", "x", "]"}]}], ",", 
    RowBox[{"{", 
     RowBox[{"x", ",", "0", ",", "l"}], "}"}]}], 
   "]"}]}], "\[IndentingNewLine]", 
 RowBox[{"C12", " ", "=", " ", 
  RowBox[{"Integrate", "[", 
   RowBox[{
    RowBox[{
     RowBox[{"phi1", "[", "x", "]"}], "*", 
     RowBox[{
      RowBox[{"H2", "''"}], "[", "x", "]"}]}], ",", 
    RowBox[{"{", 
     RowBox[{"x", ",", "0", ",", "l"}], "}"}]}], 
   "]"}]}], "\[IndentingNewLine]", 
 RowBox[{"C13", " ", "=", " ", 
  RowBox[{"Integrate", "[", 
   RowBox[{
    RowBox[{
     RowBox[{"phi1", "[", "x", "]"}], "*", 
     RowBox[{
      RowBox[{"H3", "''"}], "[", "x", "]"}]}], ",", 
    RowBox[{"{", 
     RowBox[{"x", ",", "0", ",", "l"}], "}"}]}], 
   "]"}]}], "\[IndentingNewLine]", 
 RowBox[{"C14", " ", "=", " ", 
  RowBox[{"Integrate", "[", 
   RowBox[{
    RowBox[{
     RowBox[{"phi1", "[", "x", "]"}], "*", 
     RowBox[{
      RowBox[{"H4", "''"}], "[", "x", "]"}]}], ",", 
    RowBox[{"{", 
     RowBox[{"x", ",", "0", ",", "l"}], "}"}]}], 
   "]"}]}], "\[IndentingNewLine]", 
 RowBox[{"C21", " ", "=", " ", 
  RowBox[{"Integrate", "[", 
   RowBox[{
    RowBox[{
     RowBox[{"phi2", "[", "x", "]"}], "*", 
     RowBox[{
      RowBox[{"H1", "''"}], "[", "x", "]"}]}], ",", 
    RowBox[{"{", 
     RowBox[{"x", ",", "0", ",", "l"}], "}"}]}], 
   "]"}]}], "\[IndentingNewLine]", 
 RowBox[{"C22", " ", "=", " ", 
  RowBox[{"Integrate", "[", 
   RowBox[{
    RowBox[{
     RowBox[{"phi2", "[", "x", "]"}], "*", 
     RowBox[{
      RowBox[{"H2", "''"}], "[", "x", "]"}]}], ",", 
    RowBox[{"{", 
     RowBox[{"x", ",", "0", ",", "l"}], "}"}]}], 
   "]"}]}], "\[IndentingNewLine]", 
 RowBox[{"C23", " ", "=", " ", 
  RowBox[{"Integrate", "[", 
   RowBox[{
    RowBox[{
     RowBox[{"phi2", "[", "x", "]"}], "*", 
     RowBox[{
      RowBox[{"H3", "''"}], "[", "x", "]"}]}], ",", 
    RowBox[{"{", 
     RowBox[{"x", ",", "0", ",", "l"}], "}"}]}], 
   "]"}]}], "\[IndentingNewLine]", 
 RowBox[{"C24", " ", "=", " ", 
  RowBox[{"Integrate", "[", 
   RowBox[{
    RowBox[{
     RowBox[{"phi2", "[", "x", "]"}], "*", 
     RowBox[{
      RowBox[{"H4", "''"}], "[", "x", "]"}]}], ",", 
    RowBox[{"{", 
     RowBox[{"x", ",", "0", ",", "l"}], "}"}]}], 
   "]"}]}], "\[IndentingNewLine]"}], "Input",
 CellChangeTimes->{{3.537106285738415*^9, 3.537106396522278*^9}, {
  3.537111387222766*^9, 3.537111396338312*^9}, {3.5371114374807043`*^9, 
  3.537111443202018*^9}}],

Cell[BoxData[
 RowBox[{"-", 
  FractionBox["1", "l"]}]], "Output",
 CellChangeTimes->{
  3.537106322838333*^9, 3.537111313335195*^9, {3.537111401325001*^9, 
   3.537111417491973*^9}, 3.537111448411471*^9}],

Cell[BoxData[
 RowBox[{"-", "1"}]], "Output",
 CellChangeTimes->{
  3.537106322838333*^9, 3.537111313335195*^9, {3.537111401325001*^9, 
   3.537111417491973*^9}, 3.537111448431416*^9}],

Cell[BoxData[
 FractionBox["1", "l"]], "Output",
 CellChangeTimes->{
  3.537106322838333*^9, 3.537111313335195*^9, {3.537111401325001*^9, 
   3.537111417491973*^9}, 3.537111448433526*^9}],

Cell[BoxData["0"], "Output",
 CellChangeTimes->{
  3.537106322838333*^9, 3.537111313335195*^9, {3.537111401325001*^9, 
   3.537111417491973*^9}, 3.537111448434236*^9}],

Cell[BoxData[
 FractionBox["1", "l"]], "Output",
 CellChangeTimes->{
  3.537106322838333*^9, 3.537111313335195*^9, {3.537111401325001*^9, 
   3.537111417491973*^9}, 3.537111448490192*^9}],

Cell[BoxData["0"], "Output",
 CellChangeTimes->{
  3.537106322838333*^9, 3.537111313335195*^9, {3.537111401325001*^9, 
   3.537111417491973*^9}, 3.537111448490838*^9}],

Cell[BoxData[
 RowBox[{"-", 
  FractionBox["1", "l"]}]], "Output",
 CellChangeTimes->{
  3.537106322838333*^9, 3.537111313335195*^9, {3.537111401325001*^9, 
   3.537111417491973*^9}, 3.537111448491496*^9}],

Cell[BoxData["1"], "Output",
 CellChangeTimes->{
  3.537106322838333*^9, 3.537111313335195*^9, {3.537111401325001*^9, 
   3.537111417491973*^9}, 3.53711144849227*^9}]
}, Open  ]],

Cell[BoxData[
 RowBox[{"KS11", " ", "=", " ", 
  RowBox[{"Integrate", "[", 
   RowBox[{
    RowBox[{
     RowBox[{"phi1", "[", "x", "]"}], "*", 
     RowBox[{"phi1", "[", "x", "]"}]}], ",", 
    RowBox[{"{", 
     RowBox[{"x", ",", "O", ",", "L"}], "}"}]}], "]"}]}]], "Input",
 CellChangeTimes->{{3.5371113210428457`*^9, 3.537111346734078*^9}}]
},
WindowSize->{1920, 1056},
WindowMargins->{{0, Automatic}, {Automatic, 0}},
FrontEndVersion->"8.0 for Linux x86 (64-bit) (February 23, 2011)",
StyleDefinitions->"Default.nb"
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[557, 20, 1490, 41, 30, "Input"],
Cell[CellGroupData[{
Cell[2072, 65, 267, 4, 30, "Input"],
Cell[2342, 71, 454, 9, 30, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[2833, 85, 399, 10, 30, "Input"],
Cell[3235, 97, 576, 11, 30, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[3848, 113, 120, 2, 30, "Input"],
Cell[3971, 117, 305, 6, 30, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[4313, 128, 138, 3, 30, "Input"],
Cell[4454, 133, 252, 4, 30, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[4743, 142, 141, 3, 30, "Input"],
Cell[4887, 147, 235, 3, 30, "Output"]
}, Open  ]],
Cell[5137, 153, 2654, 83, 107, "Input"],
Cell[CellGroupData[{
Cell[7816, 240, 291, 5, 88, "Input"],
Cell[8110, 247, 720, 22, 53, "Output"],
Cell[8833, 271, 650, 19, 53, "Output"],
Cell[9486, 292, 737, 23, 53, "Output"],
Cell[10226, 317, 602, 18, 53, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[10865, 340, 417, 10, 88, "Input"],
Cell[11285, 352, 623, 18, 51, "Output"],
Cell[11911, 372, 566, 16, 51, "Output"],
Cell[12480, 390, 647, 19, 51, "Output"],
Cell[13130, 411, 566, 16, 51, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[13733, 432, 120, 2, 30, "Input"],
Cell[13856, 436, 286, 5, 30, "Output"]
}, Open  ]],
Cell[14157, 444, 1026, 28, 88, "Input"],
Cell[CellGroupData[{
Cell[15208, 476, 1094, 32, 126, "Input"],
Cell[16305, 510, 381, 6, 30, "Output"],
Cell[16689, 518, 452, 9, 30, "Output"],
Cell[17144, 529, 364, 5, 30, "Output"],
Cell[17511, 536, 433, 8, 30, "Output"]
}, Open  ]],
Cell[17959, 547, 379, 9, 50, "Input"],
Cell[CellGroupData[{
Cell[18363, 560, 183, 3, 50, "Input"],
Cell[18549, 565, 188, 4, 45, "Output"],
Cell[18740, 571, 165, 3, 45, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[18942, 579, 193, 3, 50, "Input"],
Cell[19138, 584, 72, 1, 30, "Output"],
Cell[19213, 587, 70, 1, 30, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[19320, 593, 2598, 83, 183, "Input"],
Cell[21921, 678, 205, 5, 47, "Output"],
Cell[22129, 685, 184, 4, 30, "Output"],
Cell[22316, 691, 187, 4, 47, "Output"],
Cell[22506, 697, 167, 3, 30, "Output"],
Cell[22676, 702, 187, 4, 47, "Output"],
Cell[22866, 708, 167, 3, 30, "Output"],
Cell[23036, 713, 205, 5, 47, "Output"],
Cell[23244, 720, 166, 3, 30, "Output"]
}, Open  ]],
Cell[23425, 726, 344, 9, 30, "Input"]
}
]
*)

(* End of internal cache information *)
